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Abstract 

Wc study the kinematics of the Galactic thin and thick disk populations us- 
ing stars from the RAVE survey's second data release together with distance 
estimates from Breddels et al. (2009). The velocity distribution exhibits the 
expected moving groups present in the solar neighborhood. Wc separate thick 
and thin disk stars by applying the X (stellar-population) criterion of Schuster 
et al. (1993), which takes into account both kinematic and metallicity informa- 
tion. For 1906 thin disk and 110 thick disk stars classified in this way, we find 
a vertical velocity dispersion, mean rotational velocity and mean orbital eccen- 
tricity of (fjw, (V*), (e))thm = (18±0.3 km s'S 223±0.4 km s'S 0.07±0.07) and 
(o-w, (V$), (e))thick = (35±2 km s-\ 163±2 km s-\ 0.31±0.16), respectively 
From the radial Jeans equation, we derive a thick disk scale length in the range 
1.5 — 2.2 kpc, whose greatest uncertainty lies in the adopted form of the under- 
lying potential. The shape of the orbital eccentricity distribution indicates that 
the thick disk stars in our sample most likely formed in situ with minor gas-rich 
mergers and/or radial migration being the most likely cause for their orbits. 
Wc further obtain mean metal abundances of ([M/H])thin = +0.03 ± 0.17, and 
([M/H])thick = —0.51 ±0.23, in good agreement with previous estimates. We es- 
timate a radial metallicity gradient in the thin disk of -0.07 dex kpc~^, which is 
larger than predicted by chemical evolution models where the disk grows inside- 
out from infalling gas. It is, however, consistent with models where significant 
migration of stars shapes the chemical signature of the disk, implying that radial 
migration might play at least part of a role in the thick disk's formation. 
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1. Introduction 



The Radial Velocity Experiment (RAVE) is a spectroscopic survey to mea- 
sure radial velocities and stellar atmospheric parameters (Toff, logy, [M/H]) 
of up to one million stars using the Six Degree Field multi object spectro- 
^raph on the 1.2 m UK Schmidt Telescope o f the Anglo-Australian Observatory 



Steinmetz et all l2008HZwitter et~all . l2008t ). For most of the stars, RAVE also 



contains 2MASS JHK photometry, USNO-B proper motions (/i^, /i^), and es- 
tima tes of [a/Fe], althou gh the latter possess great uncertainties (typically 0.15 



dex, IZwitter et all |2008[) 



Using empirical relation s between Mk and [ J — K) of main-sequence and 
giant stars from Hipparcos, IVeltz et al. I (l2008l) have studied the kinematics of 
the thin and thick disks based on G and K type stars in RAVE. For this kind of 
work, distance s are o f great importance in order to derive 3D space velocities. 



Klement et al.l (|2008l ) took a similar approach and derived absolute magnitude 
relations between Mv^ and {Vt — H), where Vr denotes the magnitude in the 
Tycho-2 V band. They assumed that the majority of RAVE stars are main- 
sequence stars with solar metallicities, w hich has later been shown to be an 
invalid assumption ( Seabroke et al. . 2008fl . The RAVE survey contains a large 
fract ion of giants, for which photometric parallaxes can not be estimated. There- 
fore, iBreddels et al] (|2010l) published an important method to derive distances 
from the (J — K) color and the astrophysical parameters by finding the closest 
match of each star to a set of theoretical isochrones. Due to the limitations of 
their theoretical isochrone grids, the distances provided by Breddels et al. are 
only reliable for stars with log 5 > 3. 

In this paper we use a sample of ^ 4000 putative main-sequence stars (log g > 
3) from the second RAVE data release and classify the Galact ic disk popula- 
tions a ccording to the X stellar population parameter defined bv lSchuster et al 



([1993!). The X parameter utilizes both kinematic and metallicity information. 
This paper is organized as follows: § [2] describes the data, distances and selec- 
tion criteria. Our procedure to calculate the orbital parameters of our sample is 
given in § [3l the thin/thick disk classification and their kinematic and chemical 
properties are presented in § 21 our conclusions are summarized in § [S] 



2. Data, distances and selection criteria 

The second RAVE data release (DR2) allows us to derive relations between 
kinematics and elemental abundances in the Galactic disk components. We 
use a sample of 9696 putative main sequence stars from DR2, which have been 
separated from subgiants and giants by fitting Gaussian mixtur e models to the 



distri bution of spectroscopic log g estimates (for more details, see lKlement et al 



2010h . The logg values of these stars range from 3.21 < logg < 5.0 for stars 



with {J - K) > 0.5 (roughly K and M stars) and from 3.09 < logg < 5.0 fo r 
the bluer stars. We adopt the distance estimates from Breddels et al. ( 2010l ). 



because they show no evidence for clear systematic errors in these ranges of 
log (7 values. 
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We prepare our sample by putting constraints on the errors of distance, 
ra dial velocity and each of the proper m otion components. B ased on Fig. 1 
of Steinmetz et al. I (|2008l) and the work of iBreddels et alj I (l2010l) . we choose to 
remove all stars with radial velocity and proper motion component errors greater 
than 5 km and 6 mas yr~^, respectively. We further remove all stars with 
relative distance errors larger than 40%. After these cuts, there remain 4027 
putative main sequence stars for studying the kinematic and chemical differences 
between the thin and thick disks. Our stars span a range of —1.53 < [M/H] < 
+0.49 in metal abundances; the [a/Fe] values are limited to the range [0, +0.40] 
and have large recovery errors of typically 0.15 dex, because the theoretical 
template spectra used t o derive [a/F e ] are limited to only two grid points at 
[a/Fe]=0.0 and +0.40 (jZwitter et all . boOSi) . Although RAVE is not able to 
measure the a abundances of individual stars accurately, the assignment of 
[a/Fe] is not random, and it might be interesting to look at average trends of 
[a/Fe] for the different stellar populations. The galactic coordinates of our 4027 
sample stars are in the ranges of 26° < |6| < 88° and 3° < Z < 355°. 



3. Galactic Space Velocities and Orbital Parameters 



The Galactic spac e velocity components ( [/, V, W ) were computed by apply- 
ing the algorithms of Johnson fc Sodcrblom, (jl987t ) to the basic observables of 
our 4027 RAVE main-sequence stars: celestial coordinates (a, S), proper motion 
components (^q, fis), the ra dial velocity (Vj-ad) and the pa rallax tt. The trans- 
formation matri ces given in Johnson fc SoderblomI (|1987l ) have been updated 



to epoch J2000 ( Murray . 19891 ). We adopt a right-handed coordinate system 
with the (x, y, z)-axes pointing towards the Galactic center {I = 0°,5 = 0°), 
the direction of Galactic rotation {I = 90°, 5 — 0°) and the North Galactic 
Pole {b = 90°), respectively. {U,V,W) are the corresponding Cartesian com- 
ponents of a star's velocity vector with respect to the Local Standard of Rest 
(LSR). Thereby, we corrected for a solar motion with respect to th e LSR of 
(J7, V, W)q = (+7.5, +13.5, +6.8) km s'^ (|Francis fc Andersonl . [20091 ) . 

We have computed the peri- and apogalactic distances (i?niin, -Rmax) of our 
sample stars by integrating each star' s orbit for 3 Gyr in a Gala ctic potential 
consis ting of a Miyamo to-Nagai disk ( Mivamoto fc Nagal 1 19751) , a Hernquist 
bulge teernq uis^. 1990() and a logarithmic spherically symmetric halo potential 
(e.g. Johnston et al I1l99a): 



$disk(r') 



GMdisk 



\J x"^ + y'^ + {a + y/z^+ h^f 



*bulgo(r) : 

/ -r c 

$haio(r) =0.5«2 In(r2+d2) 



(1) 



Here, r — -y/a;^ + + z^, and the potential parameters have been chosen in 
order to provide a nearly flat rotation curve with circular velocity Vc ~ 220 km 
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s at the sun (see also Section lOl below): Afdisk = 1-0 x IO^^Mq, Mbuigc = 
3.4 X 101°Mq, a = 6.5 kpc, b = 0.26 kpc, c = 0.7 kpc, d = 12.0 kpc and 
vq — 180.0 km s~^. For the integration, we used a leap-frog algorithm and 
100,000 time steps, which lead to an energy conservation better than 10~® 
for most stars. Each star's {U, V, W) velocity components have been trans- 
formed into a Gala ctocentric restframe by adding a LSR velocity of 220kms~^ 
( Gunn et al. . 19791) to V. The orbital eccentricity is given by the relation 
e = (-Rmax - -Rmin)/(-Rmax + -Rmin)- Here, i?niax and i?min are apogalactic and 
perigalactic distances, respectively. For each star, we also compute its mean 
Galactocentric distance, or mean orbital radius, Rm, as the mean of apogalac- 
tic and perigalactic distances. Contrary to i?min and i^maxi Rm, remains fairly 



robust against isotropic d iffusion of stellar orb its ( Grenon . 1987t ). which is also 
referred to as "blurring" (jSchonrich fc Binnevl . 120091 ). This becomes important 
when investigating the evidence for an intrinsic radial metallicity gradient in 
the thin disk. 

The distribution of the Galactic heliocentric velocities in the ([/, V) plane is 
shown as a Bottlinger diagram in the upper panel of Fig. [1] The stars appear 
not to be smoothly distributed as expected from the presence of moving groups. 
In order to enhance any concentration of stars, we use the wavelet transform 
analysis, which results in a visual representation of the stellar density in the 
Bottlinger diagram. We construct a quadratic grid with pixels of 1 km s~^ 
width on each side and compute the wavelet coefficients at each grid point 
through 



'w{x,y) ^ / / dx' dy''^{x — x' ,y — y')y. 



N 



^5{x' - Xi)5{y' - yi) 

i=l 

N 

= ^■^{x-x^,y-y^), 



(2) 



where N — 4027 is the number of stars in our sample. For the analyzing wavelet 
^P, we use the so-called Mexican hat, which is the negative normalized second 
derivative of a bivariate Gaussian: 



^{x,y) = (2- 



y 



) 



(3) 



This type of analyzing wavelet has been used often in order to enhance over- 
densities in velocity diagrams (Skulian^etiD, Il999t lArifvanto fc Fuchd . 120061: 
Klement et al.l . l2008tlZhao et al., . ,2009,) . It is normalized such that the total area 
under its curve is zero, so that any overdense bins will have a positive wavelet 
coefficient. The wavelet coefficient in such a bin will have a maximum value if 
the scale parameter a is equal to the half- width of the overdensity, assuming it's 
'bump' is of Gaussian shape. We expect the extend of t he overdensities as soci- 
ated with moving groups to lie in the range 3-6 km s~^ ( Zhao et al. . 20091 ) and 
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Figure 1: Upper panel: The Bottlinger diagram for our sample of 4027 RAVE dwarfs. The 
displayed velocity components refer to the LSR. Bottom panel; Contours of the wavelet 
transform of the distribution of our data in {U, V) space. The overdensities correspond to 
well-known moving groups in the solar neighborhood. The color bar shows the values of the 
wavelet coefficients that follow from eq. I2I 1. 
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thus to be dominated by our typical velocity uncertainties that are somewhat 
larger. For the scale parameter cr, we therefore chose 7 km s~^, which is of the 
order of the median U and V uncertainties. The median is more representative 
for the typical velocity uncertainty than the mean because it is more robust 
against any extreme outliers. From the bottom panel of Figure [U we find the 



main moving groups expected in a solar neighborhood sample (see IZhao et "al 



2OO9I and references therein). Note that the location of many moving groups 



in velocity space is often described with respect to the sun, which we have 
account ed for in Figure [T] wh ere we plot velocities with respect to the LSR. 

Like iBreddels et al. ()201Cll ) , we place a cylindrical volume on the sun with a 
radius of 500 pc and a height of 600 pc. For stars in this cylinder, the velocity 
dispersions {au,av,crw) and mean velocities ((C/), {V), (W)) are estimated 
as (37±0.4, 26±0.3, 19±0.2) and (-11.9±0.6, -20.2±0.4, -8.1±0.3) kms-\ 
respectivel y. These values are almost exactly the same as the ones given in 
Table 1 of iBreddels et all (|2010l ). which is not surprising since our sample is a 
subset of the data used by these authors. 

We further investigate how the velocity dispersions change with respect to 
metallicity. This is shown in Table [T] and Figs. 2(a)-(c), where we tabulate 
the velocity dispersions for stars in several metallicity bins. While au and ay 
remain fairly constant over the sampled metallicity range, we observe a rise 
in aw with declining metallicity. This rise is displayed in Fig. 2(a): aw rises 
very weakly from (aw) ~ 15 km s~^ to 20 km s~^ over the metallicity range 
-0.50 < [M/H] < +0.50 and begins to rise more steeply for [M/H] < -0.50. In 
the range [M/H] = [-0.80, -0.50], (aw) is 26 km s^^. For [M/H] < -0.80, (aw) 
h as increased to 40 km s~^ . Our [M/H]-tTiv relation can be compared to Fig. 10 
of IWvse fc Gilmord (|l995l ). which shows a very similar behav ior for a sample of 



F/G s t ars drawn from the Glie se catalog and the sample of lEdvardsson et al 
( 1993f) . IWvse fc Gilmord (|l995l) note that between [Fe/H]= -0.40 and -0.50, 
there is a transition from 19 to 42 km s~^ in aw- Similarly, an abrupt in- 
crease in the vertica l velocity dispersion at [Fe/H] ^ —0.70 has been sh own by 



Gilmore et al. J 19891) (their F ig. 6b) for stars taken from the samples of NorrisI 
(|l987t ) and lStromgrenI (|l987l) . 

We interpret the trend observed in Figure 2(a) by considering two effects: 
the first is continuous heating of thin disk stars, e.g. through spiral waves, which 
gives the older stars higher vertical velocity dispersions. There is no guarantee, 
however, that older stars in our sample are indeed more metal-poor on average 
than the younger ones (e.g. iHolmberg et"all l2007h . In fact, the age-metallicity 
relation of thin disk stars is st ill under debate, and i t's shape for RAVE stars is 
currently under investigation ( Anguiano et al. . 20091) . The weak rise of (aw) in 
the range [M/H] = [—0.5, +0.5] could also be caused by a growing contribution 
of thick disk stars that have intrinsically higher velocity dispersions. This second 
efi^ect is most probably responsible for the sharp rise of the velocity dispersion 
at [M/H] < —0.80: our value of aw = 40 km s~^ is too large to be explainable 
by disk heating alone, even if we ass ume a strong correlatio n of stellar age and 



metallicity (see also the discussion in IHolmberg et al.U2007l Sec. 11). 



We point out that our stars have estimates of total metallicity rather than 
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Table 1: Metallicities ([M/H], ([M/H])) and (crjy , crv, ctw) velocity dispersions for our 4027 
putative main sequence stars. The covariances (^fjy, (^i/w '^vw ^"^^ listed in Cols. 7—9. 



[M/H] 
(dex) 


<[M/H]) 
(dcx) 


N 


"U 
km a-1 


"V 
km =-1 


km =-1 


km'^'r- 


k5r-- 


kl^r- 


1-1.53, 


-1.001 


-1.26 + 0.16 


5 


37 + 13 


59 + 21 


48±17 


1057+83 


-115 + 59 


-1356 + 55 


(-1.00 


-0.80] 


-0.87 + 0.06 


9 


45 + 11 


33 + 8 


33 + 8 


296 + 61 


-153 + 58 


245 + 52 


(-0.80 


-0.70] 


-0.74 + 0.02 


19 


48 + 8 


39 + 7 


29 + 5 


-424 + 55 


-160 + 53 


-145 + 45 


{-0.70 


-0.60] 


-0.63 + 0.03 


48 


36 + 4 


37 + 4 


24 + 2 


56 + 53 


- 8 + 43 


57+45 


(-0.60 


-0.50] 


-0.54 + 0.03 


144 


39 + 2 


27 + 2 


24+1 


166 + 51 


183 + 50 


150 + 40 


{-0.50 


-0.40] 


-0.44 + 0.03 


304 


35 + 1 


24+1 


20+1 


91 + 45 


15 + 41 


70 + 33 


(-0.40 


-0.30] 


-0.34 + 0.03 


494 


38 + 1 


24+1 


19+1 


161 + 48 


-29 + 42 


-24 + 42 


(-0.30 


-0.20] 


-0.24 + 0.03 


581 


38 + 1 


25 + 1 


19+1 


-49 + 44 


-35 + 42 


25 + 32 


{-0.20 


-0.10] 


-0.14 + 0.03 


605 


37+1 


24+1 


20+1 


247 + 49 


-10 + 42 


-17 + 31 


(-0.10 


-0.00] 


-0.04 + 0.03 


772 


39 + 1 


25 + 1 


20+1 


106 + 49 


12 + 44 


4+26 


{-0.00, 


+ 0.10] 


+0.05 + 0.03 


410 


34 + 1 


24+1 


18+1 


-25 + 41 


-48 + 37 


1 + 30 


( + 0.10, 


+0.20] 


+0.15 + 0.03 


309 


37 + 1 


25 + 1 


18+1 


41 + 46 


-28 + 40 


20 + 31 


( + 0.20, 


+0.30] 


+0.25 + 0.03 


222 


36 + 2 


25 + 1 


17+1 


127 + 47 


23 + 40 


-9 + 30 


( + 0.30, 


+0.40] 


+0.34 + 0.03 


81 


32 + 3 


27 + 2 


16+1 


168 + 46 


18 + 36 


1 + 31 


( + 0.40, 


+ 0.50] 


+ 0.44 + 0.03 


22 


35 + 5 


30 + 5 


15 + 2 


-104 + 44 


-90 + 36 


131 + 37 



iron abundances; the latter are usually lower by ^ 0.1 — 0.3 dex compared to 



[JVI/H] at low metallicities (jZwitter et al.l . l20Q8l Fig. 17). However, qualita- 



tively, our [M/Hj-criy r elatio n is in good agreeme nt with the earlier findings of 
IWvse fc GilmOTg (|l995h and iGilmore et all (|l989h . 

For the metallicity ranges from Table [TJ we compute two estimates for the 
mean rotational velocity in the Galactic rest frame: (Vj-pt) is derive d on the 
basis of radial velocities and distances alone ( Frenk fc White . 1980[ ). so that 
(eventually systematic) errors in the proper motions have no influence on the 
inferred value of rotation; thereby, the Galaxy is assumed to be axisymmetric. 
(V$) is based on the full set of observables including proper motions and is the 
mean rest-frame rotational velocity in a Galactocentric cylindrical coordinate 
system. In calculating both rotational velocities, we have assumed that the sun 
lies a t a distance Rq = 8.5 kpc from the Galactic center (jKerr fc Lvnden-Beli 
iisi), and that the velocity of the LSR is Vlsr = 220 km s ^. As can be 
seen from Table[2l both measures of rotational velocity, (Kot) and (Vb), remain 
almost constant with changing metallicity and are consistent with rotational 
velocities ty pical for the thin disk. The rotational velocities obtained from the 
algorithm of iFrenk fc Whit^ ( 1980l ) are s lightly smaller than the (V$) values. 
This behavior has already been noticed bv lCaroUo et al. ( 2007 ). who also found 
(Kot) generally smaller than (Vi) (see their Supplemental Table 2). 

From the first two [M/H] bins listed in Tables [T] and [21 two likely halo 
stars have been excluded from the statistics in order not to effect the rotational 
velocities. Their space velocities are given separately in Table [3] The mean 
eccentricities given in Table [2] decline to more circular orbits with decreasing 
metallicity. Simultaneously, ([a/Fe]) values drop to solar with some fluctuations. 



4. Separation of Galactic thin and thick disk populations and deriva- 
tion of their chemo-kinematical properties 

Separating thin and thick disk stars in the ([M/H],V ) plane 

We wish to separate our sample into members of the thin and thick disk. On 
average, thin disk stars tend to be younger, more metal-rich and faster rotating 
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Figure 2: uvc, cv, <7u versus [M/H] for the data in Table 1. 



Table 2: Mean metallicities, rotational velocities, eccentricities and alpha element abundances 
for our 4027 putative main sequence stars. The uncertainties in ([M/H]), (e) and ([a/Fe]) are 
the standard deviations. The uncertainties in V^ot and are the mean errors. 



[M/H] 
(dex) 


([M/H]) 
(dex) 


N 


km s ^ 


km s ^ 






([«/Fc]) 
(dox) 


[-1.53, 


-1.00] 


-1.26 + 0.16 




131 + 32 


219 + 26 





21 + 0-15 


+ 0-18 + 0-07 


{-1.00, 


-0-80] 


-0.87 + 0.06 


9 


117 + 30 


203 + 11 





21 + 0-15 


+ 0-19 + 0-13 


(-0.80, 


-0.70] 


-0.74 + 0.02 


19 


172 + 15 


206 + 9 





19 + 0.13 


+ 0.24 + 0.11 


(-0.70, 


-0.60] 


-0.63 + 0.03 


48 


182 + 9 


219 + 5 





14 + 0.14 


+ 0.25 + 0.10 


(-0.60, 


-0.50] 


-0.54 + 0.03 


144 


175 + 5 


215 + 2 





14 + 0.10 


+ 0.24 + 0.10 


(-0.50, 


-0.40] 


-0.44 + 0.03 


304 


185 + 3 


219+1 





13 + 0.08 


+ 0.21+0.08 


(-0.40, 


-0.30] 


-0.34 + 0.03 


494 


183 + 2 


217+1 





13 + 0.09 


+ 0.17 + 0.09 


(-0.30, 


-0.20] 


-0.24 + 0.03 


581 


181+2 


215+1 





13 + 0.09 


+ 0.13 + 0.08 


(-0.20, 


-0.10] 


-0.14 + 0.03 


605 


184 + 2 


214+1 





13 + 0.09 


+ 0.10 + 0.09 


(-0.10, 


-0.00] 


-0.04 + 0.03 


772 


180 + 2 


212+1 





14 + 0.09 


+ 0.07 + 0.07 


(-0.00, 


+ 0.10] 


+ 0.05 + 0.03 


410 


177 + 2 


212+1 





13 + 0.09 


+ 0.07 + 0.07 


( + 0.10, 


+ 0.20] 


+ 0.15 + 0.03 


309 


177 + 3 


209+1 





14 + 0.09 


+ 0.05 + 0.06 


(+0.20, 


+ 0.30] 


+ 0.25 + 0.03 


222 


181 + 3 


211 + 2 





14 + 0.09 


+ 0.03 + 0.05 


(+0.30, 


+0.40] 


+0.34+0.03 


81 


175+5 


208+3 





14+0.09 


+0.03+0.04 


(+0.40, 


+0.50] 


+0.44+0.03 


22 


197+12 


205+6 





16+0.10 


+0.01 + 0.05 
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Table 3: Distances, metallicity estimates and velocity components for the two stars in the 
first two metal abundance ranges of Tables [Hand El 



OBJECT-ID d [M/H] U V W (V"^) 
kpc (dcx) km km km km 



T8093-00436-1 0.11 -1.45 55.39 -309.82 -90.04 -90 

C0229387-130746 0.42 -0.99 -99.62 -146.29 -135.33 74 



than their thick-disk counterparts. Due to a lack of metalhcity estimates and/ or 
stehar ages, many studies h ave based their thin/ thick disk separation on the 
steUar kinematics only (e.g. |Ben sby et al., 2003). However, such an approach 
often requires fixing some values for the velocity dispersions of the thin and thick 
disk components a priori and might perform suboptimal due to some overlap in 
the velocity components of thin and thick disk stars. We therefore aim at using 
both kinematics and metallicities for our thin/ thick disk separation. 

We expect our sample to be dominated heavily by the thi n disk, because 



most of the stars possess solar-like metallicities (see also Sec. 4 in lBreddels et al 
2010( ). Due to this imbalance, we expect that any unsupervisecjl classification 



method might give non optimal results when we restrict ourselves to the RAVE 
data only. We have applied an expectation-maximization (EM) algorithm to 
the distribution of our data in the ([M/H],V$) plane. The EM algorithm fits 
a specified number of bivariate Gaussian components to the data by iteratively 
maximizing the likelihood of the data given the Gaussian model parameters. 
The EM algorithm is stopped after 100 iterations, which is sufficient to obtain 
robust parameter estimations. To prevent convergence to local maxima, we 
have first initialized the centroids of the Gaussians using a k-nearest neighbor 
algorithm with 100 iterations before each run. We checked that different initial 
guesses, that were partly motivated by the actual distribution of our data in 
the ([M/H],V$) plane, converged to the same centroids. We found that three 
components give a higher likelihood than one, two or four components. By 
connecting one of the Gaussians (say, number 1) to the thick disk component, 
we then classify star i as belonging to the thick disk if the probability given by 



^(thick I = ^3 "^::;^r^i , (4) 



is greater than 0.5. Here, 



is the "feature vector" of star i containing its metallicity and rotational velocity 
component, tt^ denotes the weight of the fcth Gaussian (with X]fc=i ""fc = 1): 
and N{xi I /.ifc, Sfc) is given as 

A/'(a;,, = —A^^ exp(-i(a;i - /Xfc)'^Sfc ^(a;» - /^fc)) . (6) 



"Unsupervised" meaning without having a labeled set of training data. 
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fik and Sfe denote the mean and (2 x 2) covariance matrix of the fcth Gaussian 
component, respectively. 

The upper panel of Figure 3 shows the distribution of our data in the 
([M/H],V$) plane; the lower panel shows contours of the three Gaussian com- 
ponents that we fit to these data. The solid red and blue contours correspond 
to 50% of the individual maximum probability density of each Gaussian, while 
the dashed ones indicate various fractions of the composite probability density 
function. The blue Gaussian is the one we attribute to the thick disk, and the 
stars classified as thick disk stars according to equation are indicated as 
blue crosses in the upper panel of Figure 3. For comp arison, we also show th ree 
straight diagonal lines that have been introduced bv ISchuster et al. (|l993l) to 



separate the thin disk, thick disk and halo stars. Each line is a linear combina- 
tion of [M/H] and V$ given by 

19.08 [M/H] + 6.63 + X 
= 0J305 ' 



where X is the stellar population parameter defined bv ISchuster et al" (1993) 



as the separation from the line X — 0, which itself is defined to run through 
th e two points ([M/H],V j,) = (-0.3,0 km s^^) and (-1.5,175 km s^^). According 
to Schuster et al. (Il993l) . the line labeled "0" divides the halo from the thick 



disk population, and stars become more disk- like with decreasing X. For our 
RAVE sample, we will classify stars as thick disk members if their location in 
the ([M/H],Kj,) plane satisfies -22 < X < --50 We see that most of the stars 
classified by our mixture model as thick disk stars are contained within this 
thick disk boundary, but many stars classified as thin disk members too. On 
the other hand, some stars are clearly misclassified by our mixture model, in 
particular three low metallicity stars with V$ < 100 km s^^ and the one star 
at ([M/H],V$)=~(0,300 km s~^). Assuming the classification based on the X 
criterion as "correct" for the moment, this means a low completeness and high 
contamination of the thick disk sample resulting from the mixture model. 

As stated before, the mixture model classification might be attenuated by 
the imbalance between thin and thick disk stars in our sample. We confirm this 
suspicion by showing that a mixture model classification of a more balanced 
data set yields more reliable results. For this, we choose the 1223 stars from 
Schuster et al. ( 2006[ ). becau se they contain a cons iderable fraction of thick disk 



and halo stars. Furthermore. ISchuster et al.l (|2006l ) have shown that a histogram 



of the X parameter for these stars can be fitted by three Gaussian distributions 
interpreted as thin disk, thick disk and halo (their Fig. 6). For a better compari- 
son to our RAVE data, we have converted the [Fe/H] values of the 



Schuster et al. 



data into [M/H] by using equation (21) from lZwitter et al.l (|2008f ): we also have 
converted their vel ocity components to t he value of the solar motion used in this 
study (taken from Francis fc Andersonl . i2009i) . We then classify them with the 



^Schuster et al.l l|l993l ) and ISchuster et al.1 l|2006t ) adopted slightly different boundaries, 
-21 < X < -6 in the ([Fe/H],V$) plane. 
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Figure 3; Upper panel: The distribution of our 4027 sample stars in the ([M/H],V$) plane. 
The red dots and blue crosses denote stars classified as belonging to the thin and thick disk, 
respectively, according to a Gaussian mixture model with three components (see text). The 
solid diagonal lines are lines of constant X; the values of X are denoted for each line. Bottom 
panel: Dashed lines indicate contours of the three component joint probability density function 
according to 0.1%, 0.5%, 2%, 10%, 50% and 90% of the maximum. Solid lines show contours 
at 50% of the maximum of each individual component. The red and blue contours correspond 
to the Gaussian componcnt{s) that we associate with the thin and thick disk, respectively. 
The horizontal dashed lines indicate rotation velocities of and 220 km s~^. 
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Table 4: Mean metal abundances, velocity dispersions, rotational velocities (Viot), (V^), mean 
eccentricity (e), and alpha element abundance ([o/Fe]) of thin and thick disk stars that are 
classified according to their X parameter. 

X N (|M/H]) <T„ a^r ^H,- (Vrot) (Vj.) <<!) (l"/F-ir 

— 22 < X < -5 TTO -0.51±0.23 57±4 36±2 35 ±2 113±7 163±3 0.31±0.16 0.19±0.11 

X < -33 1906 +Q.03ibQ.17 34±1 19±n.3 18±0.3 192±1 223ibQ.4 Q.07ifc0.07 0.12±0.Q7 



same EM algorithm we used for the RAVE sample, but now with four instead 
of three component^. We attribute two of these to the halo, and the other 
two to the thin and thick disk, respectively. A star gets classified as belonging 
to a particular class if its probability for that class exceeds the probability for 
the other two classes. Figure 4 shows the outcome of this classification which 
can be compared to a separation based on the X parameter. Red, blue and 
green crosses in the upper panel correspond to stars classified as thin disk, thick 
disk and halo members, re spectively. In the lower panel, we compare the data 
from lSchuster et ah ( 20061) (gray crosses) to our RAVE data (black dots) in the 



([M/H],V$) pl ane. For the fo r mer, typical [M/H] uncertainties at [M/H]=-0.5 
are ±1.4 dex ( Schuster et al. while for the atter we adopt ±0.2 dex 



( Zwitter et al 



2008[ ). The larger error bars of the RAVE data could be respon- 



sible for the larger spread in the [M/H] direction observed in figure 4. In the 
V$ direction, the uncertainties of both data sets are co mparably smal l (±6.0 



km s~^ for our RAVE data, ±6.6 km s~^ for the data of Schuster et al.h . Con- 
cerning the classification of stars into the Galactic components, we see that the 
boundary between thin and thick disk resulting from the mixture model can 
be approximated by the diagonal line corresponding to X = —22. Also, the 
line X = —5 agrees well with the mixture model boundary between thick disk 
and halo in the metallicity range —0.8 < [M/H] < 0.4, while there is some 
scattering of a few mixture model "halo" stars across this boundary at lower 
and higher metallicities (and vice versa, we find some "thick disk" stars in a 
region ascribed to the halo by the X parameter). Accounting for the expected 
differences between the linear decision boundaries from the X criterion and the 
more complex and curved ones from the mixture model, we conclude that both 
classification methods agree well. We use this fact together with the simi lar 
distribution of our RAVE stars and the thin disk stars from ISchuster et al.l in 
the ([M/H],V$) plane to justify the use of the X parameter for classifying our 
sample. 

4-2. Deriving the chemo-kinematical properties oj thin and thick disk 
4-2.1. Rotational velocities 

For the 110 stars classified as members of the thick disk in the previous 
section, we find a mean metallicity of ([ M/H]) = —0.5 1 ± 0.2 3, which is in 
good agreement wit h the values —0.50 of ISchuster et al. ( 19931 ) and —0.55 of 



Karatas et al. ( 2005f ). The velocity dispersions are (au, (tv,(^w) = (57±4, 36±2, 



^The likelihood for a four component fit is higher than for only three components. 
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Figure 4: Upper panel; The distribution of 1223 stars taken from [Schuster et al. 1 lll993l l 
the ([M/HJjV"^) plane. Red, blue and green crosses denote stars classified as belonging to the 
thin disk, thick disk and halo, respectively. This classification is based on a Gaussian mixture 
model with four components. The solid diagonal lines are lines of constant X with the value of 
X denoted for each line. Bottom panel: Dashed lines indicate contours of the four-component 
joint probability density function according to 0.1%, 0.5%, 1%, 2% and 10% of the maximum. 
Solid lines show contours at 50% of the maximum of each individual c lass f two Gaussians for 
the halo). The grey crosses are the 1223 stars from [Schuster et al.l ||1993| ). while the black 
dots give the location of our RAVE stars. Typical error bars of both data at [M/H]=-0.5 are 
given in the bottom left. Note the good agreement between our mixture model and the X 
criterion in separating the thin and thick disk. 
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Table 5: Literature values for the velocity ellipsoid of the thick disk. All values are given in 



units of km 


s-i. 










(JV 


<JW 


(^*> 


reference 


50±3 


56±3 


34±2 


190±5 


Chiba & Beers (2000) 


63±6 


39±4 


39±4 


169±5 


Soubiran et al. (2003) 






32±5 


154±6 


Karatas et al. (20051 


74±11 


50±7 


38±7 


178±8 


Vallenari et al. (2006) 


53±2 


51±1 


35±1 


182±2 


CaroUo et al. (2010) 


57±4 


36±2 


35±2 


163±3 


this study 



35±2) kms~^. The vertical velocity dispersion aw = (35±2) kms"-'^ of the thick 
disk agrees well with recent estimates that we summarize in T able El It als o 
agrees well w i th old er lit erature values o f 30-3 9 km s^ given bv lNorri!j(|l987h . 
Carnev et al.1 (jlQS Q') and' Croswell et alJ (|l99lh . 



For the thin disk, we find the mean metallicity to be consistent with solar: 
[M/H]= +0.03 ±0.17. The velocity dispersions are (tify, try, ctvk) = (34±1, 
19±0.3, 18±0.3) km s"^. Our estimate for the vertical velocity dispersion, 
aw = (18ib0. 3) km s~^, is in v ery g ood agreement with the value of (18±1) 
km s~^ from | Nordstrom et al. (2004) and the value of (20±1) km s^^ from 



Soubiran et al.l (|2003D 



The mean rotational velocity (V$), that we obtain for the thick disk, is 

(163±3) km s~^ and agrees well with literature values of ( 160+30) km s~^ 

(Norris','l986') and (157±4) kms'^ (jAUende Prieto et al.l .|200i). although lKaratas et al 
(20I).5) found a slightly lower value of (154±6) km . However, as is evident 
from Table m the (K-ot) value of 113±7 km s~^, which has been calculated 
without using proper motions, is smaller than that of (Vi) and does not agree 
with other estimates from the literature. This indicates that utilizing the full 
available information is superior to replacing measured informations (proper 
motions) through assumptions (an axisymmetric potential). 

4-2.2. The eccentricity distribution 

Thick disk stars in our sample display a mean eccentricity of (e) = 0.31±0.16 
(Fig. 5), while for thin disk stars we find (e) = 0.07±0.07. The latter result is 
fully expected, since it implies almost circular orbits for thin disk stars. The 
mean eccentricity we find for the thick disk supports the view that accretion as a 
major mechanism in thick disk formation can be ruled out, at least for our stars 
at moderate heights above and below the plane. The reason is that accreted 
stars would broaden and shift the eccentricity distribution towards higher values 



([Sales et al.l . I2009i ). To gain further insights, we have overplotted i n Fig. 5 the 



xedictions of two other thick disk formation models taken from Sales et al 



2009"), namely radial migration and minor mergers. These are eccentricity 
distributions predicted for disk stars in the height range 1 kpc ^ ^ 3 kpc 
in order to minimize any contamination from the t hin disk in these models 
(assuming the thick disk scale height to be ^ 1 kpc, see lSales et al. 1,12009", Fig. 3). 
Our stars are placed much closer to the plane (z < 500 pc), but through the X- 
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criterion we neverthele ss are able t o pick out stars from the thick disk. Although 
the potential used by I Sales et aL (l2009l) to compute eccentricities is different 
from ours, both potentials result in a simil ar rotation curve nea r the sun and 
they should yield similar values of e (see also Dierickx et al. . 2010l) . We therefore 



can at least make a qualitative comp arison of our eccentricity distribution with 
the ones shown in Sales et al. I Hooi). Such a comparison suggests that minor 
mergers and/ or radial migration might be the main mechanisms responsible 
for placing our t hick disk stars at their orbits, although the radial migration 
scenario used by ISales et al. I (l2009l) predicts a narrow peak at low eccentricity 
which is not observed. 

We have further divided the stars formed in minor mergers into in situ and 
accreted st ars and overplot th e distribution of the in situ stars in Figure 5 (again 



taken from Sales et al 



The shape of this distribution fits even better 
to our observed one as it substantially decreases the number of stars in the 
high eccentricity tail (which result from accretion) . Therefore, we can conclude 
that at least the fraction of thick disk stars close to the plane have formed in 
situ. Their e-distribution suggests that a majority of them might have formed 
in situ from g as that has been a ccreted during a period of minor mergers 



at high redshift PBrook et al.l . 120041 ) . Two recent studies specifically aimed at 



inferring the origin of the thick disk from its eccentricity distr i bution have come 
to the same conclusions ( Dierickx et al. . 2010t Wilson et al.l . I2OIOI) . We note, 
however, that care must be taken with our conclusions because they are based 
on a comparison of the e distribution in different heights ab ove /below the plane. 
Furthermore, as also pointed out by lPierickx et aL ( 2010[ ). we only compare to 
one particular simulation f or each of the thic k disk formation scenarios, namely 
the one shown in Fig. 3 of Sales et al] ( 2009t ). 



4-2.3. [a/Fe] abundances of the thick disk stars 

Further clues about the formation of the thick disk could possibly be gained 
by investigating a element abundances; the a elements provide important clues 
about the history (and in particular the duration) of star formation based on 
the different explosion timescales and enrichment products of supernovae II 
and la. Furthermore, the scatter in [a/Fe] among thick disk stars gives clues 
about the mixing state of the inter-stellar medium (ISM) from which these stars 
formed. However, as already stated above, RAVE's measurements of [a/Fe] are 
not reliable for individual stars. We nevertheless have estimated mean values for 
both disk components, because star-to-star variations in [a/Fe] are not random, 
so that we at least expect to observe a correct trend. 

As can be seen from Tabled the thick disk stars have ([a/Fe]) = 0.19±0.11, 
slightly higher than thin disk stars with ([a/Fe]) — 0.12±0.07. Despite the large 
uncertainties, the mean [a/Fe] v alues of thick and thin di sks are quite s imilar 
to the values given in Fig. 1.4 of lNissenI (|2004t ). Recentlv. iRuchti et al.l (|2010l ) 
derived separate abundances of four a elements (Mg, Si, Ca, Ti) and iron for 
243 metal-poor ([Fe/H]< -0.5) stars from RAVE of which ~ 40% have been 
classified as thin or thick disk stars based on their positions and velocities. These 
authors found [a/Fe] to be enhanced in thick disk stars with low scatter over 
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Figure 5: The eccentricity distribution of our 110 thick disk stars (red filled histogram). 
Overplotted are the predictions from two thick disk formation scenarios: radial migration 
(dashed line) and minor gas-rich mergers (solid line). The subset of stars that formed during 
minor mergers in situ is shown as the dotted line. 
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a wide range of [Fe/H], although no mean values are given. This would imply 
rapid star formation on a timescale < 1 Gyr and a high degree of mixing of the 
gas from which these stars formed; this in principal is consistent with the in 
situ formation of the thick disk from accreted gas, although such a mechanism 
would predict further stars that have been directly accreted at la ter times and 
hence posess lower [a/Fe] values, but have not been observed by iRuchti et al 



( 201Ci ). We note already that our sample contains such "missing" stars at low 



metallicities, but due to their small number and the great uncert ainties on our 



a/Fe] values, we make no further comparison to the sample of iRuchti et al 



(201 



4.3. Derivation of the scale length of thick disk 

From our estimates for the asymmetric drift, V\ag, and the components of 
the velocity ellipsoid, we are able to estimate the scale length of the thick disk 
from the radial component of the steady-state Jeans equation: 

Here, vr, and Vz denote Galactocentric cylindrical velocity components, 
p{R, z) is the volume density of the thick disk stars, i.e. the tracer popula- 
tion, and $tot(-R, z) the total gravitational potential of the Galaxy. We assume 
that our thick disk stars represent a relaxed population of stars in equilibrium 
with the total gravitational potential ^totl^, -z), which is given through equa- 
tion ([1]). Close to the midplane at the sun's position, this potential is domi- 
nated by the thin dis which justifies our choice of a single disk component 



(see also Sect. 4.1 in Girard et al. . 20061 ). On top of that, we are free to choose 



an independent description for the thick disk density, for which we assume an 
exponential profile, because it allows an easy solution of the Jeans equation. 

We evaluate the radial Jeans equation for a steady state disk at i? = Rq 
and 2 = 0. This is justified by the small distances that our stars probe (2 < 700 
pc, \R — Rq\ < 950 pc). Thus we can replace the term R^§^ in equation ^ 
with w^, the squared circular velocity at Rq. 

Usually, one assumes that there is no net radial and vertical stellar motion 
in the solar neighb orhood. This assumption has recently been challenged by 
Siebert et al.l ( 2011 ). who found a radial velocity gradient of at least 3 km s^^ 



kpc"-'^ towards the Galactic center from a sample of 213,713 RAVE stars with 
l^l < 1 kpc. However, before this result is confirmed by more studies, we will 
here keep the conservative approach and use {vr) = and (vz) = 0, so that the 
term (vr) can be replaced by aj^ and (vrVz) by aRz- We further have ff = 
at 2; = by symmetry. With (w?) = + (v^)"^, the radial Jeans equation then 



^The density of the Miyamoto-Nagai disk at Rq = 8.5 kpc is 0.16 Mq pc 3, which is 284 
times larger than the contribution from the bulge and halo. 
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reads 



2 / \2 2 2 R9{pal) daR^ 



dR dz 
R djpal) R dan,- 



R P'^R '^R 

The cross term ^-^f^ vanishes if the principal axes of the velocity ellipsoid are 
aligned with the cylindrical coordinate system; our data indicate that this is not 
strictly the case. We measure a tilt angle between the long axis of the velocity 
ellipsoid and the sun-center line of a^z = 5.6° ± 0.4°, where a^z is given by 

2 

tan(2afl^) = ^ , (10) 

^R '^2 

and the uncertainty has been evaluated by computing 1000 Monte Carlo values 
using the velocity uncertainties for individual stars. Our measurement means 
that the principal long axis of the velocity ellipsoid points slightly above the 
Galactic center. This angle is small enough, however, to neglect it's influence 
on our evaluation close to the midplane. We therefore neglect the cross term in 
equation ([9]). 

The term — ^^U^ can be simplified by assuming that the shape of the 
thick disk's velocity ellipsoid is not varying with respect to R, an assumption 
that is reasonable given the small spatial extend of our data. This implies 
ct|j oc . If the thick disk would be self-gravitating, the squared vertical velocity 
disp ersion would be pro portional to the thick disk density. This would imply 
(e.g. iGirard et all |2006|) 



R d{pal) _ 2R 
p4 dR hn ^ ' 

Here, the density distribution of the thick disk is given by p(i?, 0) = po exp(— 
However, above we have stated that our thick disk stars are presumably dom- 
inated by the thin disk's gravitational potential; therefore it seems more rea- 
sonable to assume that the squared vertical velocity dispersion of the thick disk 
stars is proportional to the density of the thin disk. We thus have aj^ oc a'^ oc 
Pthin(-Ri 0), where /Othin is g iven analytically by the Miyamoto-Nagai disk density 
(jMivamoto fc Nagail Il975l ). It then follows that 



R d{pajf) _ R R dp, 



thin 



paj^ dR Hr pthin dR 



(12) 



Equations pT|) and (|12p represent two extreme cases of a totally non-self- 
gravitating and a totally self-gravitating thick disk, respectively. The truth 
presumably lies somewhere between these two cases (jGirard et al l l2006[ ). but 



we will first use equation (|T2t for self-consistency. We will also investigate, 
however, how our estimate of the thick disk scale length changes if we apply 
equation (fTT]) instead (i.e. if we assume a self-gravitating thick disk) . 
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We now evaluate equation ([SJ with our measurements and assumptions to 
obtain an estimate for /i^, the radial scale length of the thick disk. Setting 
o'R = ajj, o'tp = o-y and Vc — (w^) = Ving, equation (jH]) yields 



4 



1 



R dp, 



thii 



Pthir 



OR 



R 

hR 



(13) 



Vc is related to the Galactic potential through — R 



dR 



The potential ^ 



that we used to integrate the orbits yields a circular velocity at the solar radius 
i?0 = 8.5 kpc of 223.1 km s~^, not too different from our adopted LSR velocity 
of 220 km s-i. 

With the values for our thick disk velocity ellipsoid, {cru,crv, (w^)) = (57, 
36, 163) kms~^, and our adopted value for the the solar radius, Rq — 8.5 kpc, 
we estimate the scale length of the thick disk as 



Hr — 1.5 ± 0.3 kpc, 



(14) 



where the uncertainty is determined by 1000 Monte Carlo computations of Hr, 
each time adding a normally distributed random error to the observables with 
standard deviation given by their uncertainties in Table SI Using our adopted 
LSR velocity, Vlsr = 220 km s~^, instead of Vc changes this result slightly to 
1.6±0.3 kpc. 

Our estimated value seems too small compared to scale length estimates 
from the literature that typically range from 2.2-3.6 kpc ( Morrison et al.l . Il990 ; 



SoubiranL[l993[|Robin etall . 1 1996115] ha et al.l . [l99llVallenari et al.Ll200aljuric et al 



20081 : Carollo et al. . 2010 ). Our estimate could be biased by some of the assump- 
tions we made. One drawback is that the gravitational potential that dominates 
the motion of the tracer stars has to be globally well-tuned. This follows directly 
from equation ()13p , which depends mainly on the circular velocity and the mea- 
sured velocity dispersions. For example, by doubling both the scale parameter 
a and the bulge mass Afbuigo in our total potential ([T]) , and - not unrealistically 
- assuming the same observations, we would obtain similar values for the radial 
velocity and scale length (219 km and 1.3±0.2 kpc, respectively), despite 
having created an overall different galaxy with the local stellar density reduced 
by 34%. As another example, if we assume the thick disk to be self-gravitating, 
i.e. using equation (fTTI) . our result would change to hft = (2.2 ± 0.3) kpc. Al- 
though this seems a less likely assumption, iGirard et al.l (|2006l ) also found that 
their best thick disk models where those that assumed the "self-gravitating form 
of the pressure term for the thick disk" . 

We therefore conclude that our data together with the adopted Galaxy model 
and assumptions outlined above allow us to derive a rough estimate for the thick 
disk scale length in the range 1.5 — 2.2 kpc, which is at the lower end of other 
values estimated so far. 



4-4- The metal abundance gradient 

The relation between the maximal height above or below the galactic plane 
that a star reaches, Zmax, and the metal abundance for our thick disk stars is 
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Figure 6: [M/H] versus ^max for thick disk stars with — 22 <C X <C — 5. 



shown in Fig. 6. No vertical abundance gradient is evident based on the mean 
values of [M/H] and Zmax (filled dots) in Fig. 6. 

We also investigate the evidence for a radial metallicity gradient in the thin 
disk by using the mea,n orbi tal radii, i?m, that are fairly robust against stel- 
lar migration (|Grenonl . Il987l) . In Fig. 7, we show the trend of [M/H] versus 
Rm for our thin d isk st ars; the fitted radial gradient is — 0.07±0.01 dex kpc~^. 
Nordstrom et al.l (|2004l ) found a radial iron abundance gradient of —0.09 dex 
kpc~^ for stars younger than 10 Gyr. Our estimate is also consistent with 
a gradient of —0.10 dex k pc~^ predicted by the chemical evolution model of 



Schonrich fc BinnevI (|2009l ) and qualitatively with the gradient expected for 
stars close to the p lane (and hence young) from simulations carried out by 
Roskar et al.l (|2008l Fig. 2), although the latter authors stated no values. In 
both of these simulations, the thick disk of the Galaxy is explained naturally 
through stellar mig ration, leading t o a ne t outward and inward motion of stars 



with time (see also iLoebnian et al.l . 120101). In contrast , the ch emical evolution 



models of Chiappini et al.l ( 1997t ) and IChiappini et al.l (2001), where the disk 
is decomposed into radial annuli that neither exchange gas no r stars, predict a 
much lower radial gradient of the order ~ 0.01—0.03 dex kpc~^ (jChiappini et al 



19971 Table 4). In these "Two Infall" models, the thin disk slowly forms inside- 



out after the thick disk and halo from a separate gas infall epoch. A gas density 
threshold for star formation is applied, preventing the growth of abundances 
and the build-up of a steeper metallicity gradient in the outer regions of the 
disk (including the solar neighborhood), where the gas density is low. Thus 
it seems that radial migration models are better suited to account for the ob- 
served radial metallicity gradient than disk formation scenarios where stars are 
basically confined to their radial birth places. 
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Figure 7: [M/H] versus the mean Galactocentric distance Rm for thin disk stars with X < —33. 



5. Conclusions 



Based on a sample of 4027 putative main sequence stars selec ted from the 
RAVE second data release which have distance estimates taken from lBreddels et al 
( 2010 l). we successfully reproduce the mean kinematic properties of the thick 
and thin disk populations compared to previous results from the literature. 
The thick disk dominates our sample for metallicities [M/H]< —0.8 (Fig. 3), 
which explains the rapid rise in the vertical velocity dispersion occurring at this 
metallicity (Fig. 2a). On the other hand, we observe only a weak rise of aw 
in the range —0.50 < [M/H] < +0.50; this rise of aw could be explained ei- 
ther by the growing fraction of thick disk stars or heating of old thin disk stars 
or both. Contribution of the heating effect would imply that there is a trend 
in the thin disk's age-metallicity relation, namely that older stars are being 
on average more metal-poor. The literature on the age-metallicity relation of 
the t hin disk has yielded mix ed results, with so me data showing a cle ar trend 



(e.g. Edvardsson et al. . 1993 ). while others not ( Holmberg et al. . 2007[ ). and it 



seems that selecti on effects in th e different surveys are responsible for these dif- 
ferences (Holmb erg et al.L 12007 ) . A st udy of the age-metallic ity relation with 
RAVE stars is currently undertaken bv lAnguiano et al~ ( 20091 ) . and the reader 
is referred to their poster for a short review and outlook on this subject. 

A separation between thin and thick disk is achieved by combining kine- 
matics an d metal abundances in the so-called X stellar population parameter 
defined bv 'Schust er et al" (1993). We have estimated the velocity ellipsoids of 
the thin and thick disk to be {au,crv,crw) = (34, 19, 18) km s~^ and (57, 36, 
35) km s^^, respectively. The vertical velocity dispersion aw — 35km s~^ for 
the thick disk is in good agreement with recent literature values in the range 
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32-39 km (Table [SJ. It allows us to derive a thick disk scale length in the 
range of hft = 1.5 — 2.2 kpc, although this has to be considered a rough esti- 
mate, because the global form of the gravitational potential that dominates the 
motion of the thick disk stars is not constrained well. 

There are four principal formation mechanisms for the thick disk: (i) a 
puffed-up ("heated") thick disk resulting from the collision of a sate llite with a 
pre-existin g thin disk dOuinn et. al.l 1993 : Velazquez fc White . 19991) : (ii) radial 
migration (jSchonrich fc Binnev . 20091) : (iii) formation of thick disk stars out of 
accreted gas during a gas-rich merg er ( Brook et al.l. 2004 1; (iv) direct accretion 
of stars into the disk ()Abadi et al.l . l2003l ). I Sales et all (|2009() have compared 
the eccentricity distribution of thick disk stars resulting from each of these 
four processes and found that a prominent peak at low eccentricity is expected 
for the heating, migration and gas-rich merging scenarios, while the accretion 
scenario predicts a large number of stars on highly eccentri c orbits. As can b e 
seen from a comparison between our Fig. 5 and Fig. 3 of I Sales et al. I (l2009l) . 
our eccentricity distribution centered at (e) = 0.31 ± 0.16 matches a scenario in 
which gas-rich minor mergers and/or radial migration are the key drivers of thick 
di sk formation . We point out, however, that our stars and the stars displayed 
in Sales et al. I (l2009t) occupy different heights above and below the plane. It 



is expected that the number of stars on highly eccentric orbits increases with 
\z\, broadening the high-eccentricity tail of the e-distribution. This would make 
radial migration more unlikely than minor mergers because the former predicts a 
very narrow peak at low e and lacks a high-eccentricity tail. However, our radial 
metal abundance gradient of —0.07 dex kpc~ ^ is consistent with the chemical 
evolution model of Schonrich fc Binnev ( 20091 ). in which radial migration alone 
is sufficient to explain the chemo-kinematic separation between thin and thick 
disks. We caution again that the stars in our sample do not cover a sufficiently 
large part of the thick disk (in particular in Zmax) in order to draw definite 
conclusions on the formation of the thick disk as a whole. In addition, it would 
be interesting to compare other diagnostics of the different models, e.g. radial 
or vertical metallicity gradients. 

We have not detected any vertical metallicity gradient in the thick disk, 
while radial mixing would predict [ M/H] to fall with incre asing z due to a net 
upward motion of stars over time (jLoebman et al.l . 2010t ): however, to what 
extend strongly depends on the enrichment history of the Galaxy and on the 
undetermined inside-out formation (Ralph Schonrich, 2010, private communica- 
tion). In addition, our data are sparse at large z and exhibit a strong dispersion. 
This is a general problem in the literature and the main reason that no con- 
sensus on a vertical metallicity gradient exists yet. In fact, the large dispersion 
observed in all data sets so far may reflect the intrinsic complexity of Galactic 
evolution inc luding processes such as in-fall, out-flow, radial mixing and stellar 



interactions ( Cui et al. . l2003h 



Our value for the radial metallicity gradient of the thin disk, on t he other 

hand , is consistent with estimates from the Geneva-Copenh agen survey (iHolmberg et al 



2007 ) and predictions made by radial migration models ([Schonrich fc Binnevl . 
20091 ). These gradients are steeper than expected for the chemical evolution 
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Table 6: Same as Table |4] but now the stars have been classified into thin and thick disk 
members by using the calibrated metallicities. Note that [M/H] now refers to the calibrated 



metallicity. 
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([M/H\) 


"U 


"V 


"W 






( = > 








(dox) 


km 


km s^l 


km s^l 


km a-1 


km 




(dox) 


-22 < X < -5 


34 


-0.28±0.25 


72±9 


47±6 


49±6 


66±15 


132±8 


0.45±0.21 


0.12±0.12 


X < -33 


3216 


+ 0.12±0.16 


35±0.44 


20±0.25 


18±0.22 


187±1 


220±0.35 


0.12±0.07 


O.lOiO.lO 



models of IChiappini et al. l (ll997Ll200ll) . who assume a gas density threshold for 
star formation and a gradually declining inside-out growth of the Galactic thin 
disk. Thus it seems that radial migration to some extend is needed to account 
for the observed gradient. 

In this study we have used the uncalibrated metallicities given in the RAVE 
DR2 catalog. When using the calibrated ones instead, but keeping the definition 
of the X parameter, we only end up with 34 thick disk stars. This is expected, 
because the calibration of the metallicity parameters that are derived in the 
RAVE pipeline through a penalized comparison with a grid of synthet i c spec - 



tra shifts these on average to higher values (see Fig. 17 in lZwitter et al.l . 120081) . 



In our case, the mean [M/H] of the thick disk stars changes from —0.51 to —0.28 
when using calibrated [M/H] estimates (see Tabled. The calibration mostly 
affects metal-poor stars, so that our estimate for the radial metallicity gradient 
in the thin disk changes from — 0.07±0.01 dex kpc^^ to — 0.04z fc0.01 dex kpc 



which in pri ncipal is consistent with the "Two-Infall" model of IChiappini et al 



(|l997l I2OOII ). We also note a rise in the velocity dispersions and a drop in the 
rotational velocity of the thick disk from 163±2 to 132±8 km s^^, while the 
results for the thin disk remain stable within the uncertainties, ajj and ay 
generally show much less agreement between different studies in the literature 
(see Table [5]), so that we can not infer whether the values obtained with the 
calibrated [M/H] estimates are more reliable than the ones obtained with the 
uncalibrated metallicities. Inserting the values of (tc/, cry, Mag) from Table [6] 
into equation (|13|) . we obtain a radial scale length of the thick disk of (2.3±0.9) 
kpc, which is larger than our result for the uncalibrated metallicities and some- 
what more consistent with the literature values cited in Section|4l crwi however, 
is more consistent with recent literature estimates (Table [S]), which are all < 40 
km s~^, when we adopt the uncalibrated metallicities as done in this study. 
With the calibrated [M/H], the mean eccentricity of the thick disk stars rises 
from 0.31±0.16 to 0.45±0.21, which would speak even more strongly against 
a formation of the thick disk by radial migration alone, although the uncer- 
tainties are quite large due to the small number of stars. We conclude that by 
using the calibrated [M/H] values given in the RAVE catalog, we would end up 
with fewer thick disk stars, but our main results would not change significantly 
within the uncertainties, apart from a change of the thick disk velocity ellipsoid 
from (cr,7,cry,crvF) = (57±4,36±2,35±2) to (72±9,47±6,49±6) and an expected 
decrease in the radial metallicity gradient of the thin disk. 
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